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We evaluate vertex corrections to Tc on the basis of the antiferromagnetic spin-fluctuation 
model of the high-Tc superconductivity. It is found that the corrections are attractive in the 
d^2_y2 channel, and they become appreciable as we go through an intermediate-coupling regime 
of Tc — lOOK, the maximum Tc attainable in the one-loop Eliashberg calculation. 
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As a model for the high-Tc superconductivity, the spin 
fluctuation jpechanisni has been one of the most widely 
discussed. Bd' The model assumes that quasiparticle is 
coupled with antiferromagnetic spin fluctuation, repre- 
sented by a peculiaii Jow-energy expression for the mag- 
netic susceptibility.EfB' The phenomenological coupling 
constant to fit the transition temperature Tc is used to 
explain, among others. plhe_ a|i|U3|raaloU|S-,transp^ prop- 
erties consistently.Eiy'iyQ§ByBBy) On the other 
side, from a microscopic point of view, numerical studies 
based, on the fluctuation exchange (FLEX) approxima- 
tionlil-' have been carried out by many authors to esti- 
mate Tc as well as to explain_.thC|-dpaations from the 
normal Fermi liquid behavior. til't3ll3'll3) Computational 
feasibility of these strong coupling theories rests on the 
effective use of a fast Fourier transform (FFT) algorithm. 
To asses the quantitative aspect of the theories, correc- 
tions coming from higher-order terms are investigated 
for the vertex function at some flxed external momenta, 
e.g., on the basis of the spin-fluctuation model, and the 
qualitativa^ featiims , nf. the effect are emerging to some 
extent.yyMBBBES However, the total effect of the 
vertex corrections on the physical observables is yet to 
be estimated numerically. Indeed to do this is gener- 
ally formidable because of the inapplicability of the FFT 
to a required additional sum on internal frequency and 
momentum. In this paper, we manage to evaluate the 
vertex corrections to Tc, and discuss questions of conver- 
gence of the formal perturbation theory with respect to 
the coupling constant. 

To put it concretely, the following investigatiai 
based on the model of Monthoux and Pines (MP),Q' 
in which the self-energy S(p, ia;„) is determined as a self- 
consistent solution of the equations. 
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G(p,ic^„)-i = G(")(p,ic^„)-i - E(p, ia.„) + (5m, (2) 



where 



and 



G(°)(p,ic.„) = 7- 



1W„ 



x(o)' 



(3) 



Ep — ~2t(cOSPx + COS Py) — At' cos Pa; COS Py. (4) 

In eq. (P, x(q, ii^m) as a function of the Matsubara fre- 
quency I'm = 2TO7rT is inferrecLfrom the low-energy form 
of the magnetic susceptibilityP-* 



x{q,^) 



(5) 



where 



LOsfil+eiq-Qf), g = (^,7r), (6) 



for Qx > and Qy > 0. We assume 
1 r° lmx{q,uj) 



XiQ^t^n) 
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X{q,uj^0). 



(n^O) 

(n = 0) (7) 



Here it is noted that a cutoff loq has to be artificially 
introduced so as to meet the condition x('7,i^'n) -^ ^/^n 
as \i'n\ -^ oo.q' For eq. (^), the integral in eq. (7) is 
analytically evaluated; 
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2XQ^sf/7r 



k^n tan 
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(8) 

To estimate the critical temperature T^ the linearized 
gap equation is used, 

T 



^(p,iw„) = --^Yl ^(P:i^n;p',iw«') 



p' ,ri' 



x|G(p',iu;„OI^^(p',i'^n'): 



(9) 
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where $(p, ia;„) is the anomalous self-energy. The pair- 
ing potential V {p, iun] p' ,'vjJn') reads 

y(p, itJ„;p',icJ„') = V^'^\p - p' ,lUJn -ic^ri') 
+ V^''\p,iUn;p' :iUJn') +VP{p,iLOn\p' ^yUJn'), (10) 
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where 

V^^Hp-p^iuJn -iuj,,,) = g'^x{p-p'^^n -iw„')- (H) 

The second and third terms in eq. (10) originate from 
the vertex corrections that we discuss below. 

As we are concerned about the d-wave instabihty, in- 
troducing a notation 



{f{p))p = 7^ Y^{C0S{P^) - COs{py))f{p), 



we put eq. (9) into 



] K(iu!n,iuJn')<l'iiiUn'), 



(12) 



(13) 



where 



and 



<P{iujn)-- {•P{p,iu;n))p- 






(14) 
(15) 



xp) = -r(y/^'''|G(p',ic.„oi') ,■ (16) 

/ON 

Here K^ (iw„, iuJn') {i — v,c) conic from the vertex cor- 
rections. As is clear from eq. (Esl), the condition that 
the largest eigenvalue of K{iujn, '^^n') reaches unity pro- 
vides a nonzero solution ^(iti;„), thus defines T^. Below 
we look for a real solution <?(ia;„), for which the imagi- 
nary part oi K{iujn,\ujni), namely, lTaV{p,iujn;p' ,\ujn'), 
is neglected. 

As for the parameters, we assume t = 0.25eV and 
t' = -0.45i for eq. (0), and we take 



e = 2.3, XQ = 75/eV, uj^i = 14meV, 



(17) 



to describe_\;(q,a;) of YBa2Cu307 at Tc = 90K, accord- 
ing to MP.EP The chemical potential /i^"^ is fixed by 

p 

(18) 

in which we assume n = 0.75 throughout this paper. The 
shift (5/i in eq. (0) is adjusted in every iteration to assure 



p,7l 



2T 



Sn=—J2 [GiP, ic^n) - G(°) (p, iu;„) 



0, 



(19) 



p,n 



as this is easier to handle than the formally divergent 
sum, n — [2T/n) ^ ^ G{p, ia;„). For practical purposes, 
the Matsubara sums in eqs. (|l|), (9) and (19) are re- 
stricted within a finite range Iw^l, jt'ml < ujc- To avoid 
spurious temperature dependences, the cutoff is fixed at 
uJc = 6.2eV ~ 3 times the bandwidth, for which we have 
|i/„| ~ LOc for m = ±2^ at T = 90K& 
The above, except for the cutoff wq in eq- (7), are all thfi 
necessary ingredients to reproduce the results of MP.LF 
Nevertheless, we could not derive them precisely, though 
qualitative features are consistently reproduced. Let us 
discuss the point briefly. First we have to note the ujq- 
dependence of Tc as a function of the coupling constant 
g^. In Fig. 1, g^ to give Tc = 90K is shown as a func- 
tion of LUo- In this result, K'^^^ in eq. ( [l6| ) is used for 



Fig. 1. The coupling constant g^ to account for Tc = 90K as 
a function of the cutoff luq, calculated on a 16x16 lattice with 
periodic boundary conditions. 



Fig. 2. Critical temperature Tc is shown as a function of g^ 
for the cutoff luq = 0.4eV. The other parameters are given in 
the text. The results are shown for a 16x16 (circles), 32x32 
(squares), 64x64 (diamonds) and 128x128 (triangles) square 
lattice. The effect of E(p, iaj„) is not taken into account for 
the open symbols, while it is included for the closed symbols. A 
closed triangle for 128x128 is not shown. 



is not taken into account, but the effect of I](p, iwn) is. 
As we see from the figure, we cannot fix g^ just from 
Tc{g^) — 90K without the knowledge of the cutoff uiq in 
eq. (7). The strong dependence on coq in the low-energy 
region ujq < 0.3eV reflects that the transition tempera- 
ture Tc, unlike the transport properties controlled by the 
quasiparticle damping, is not determined solely from the 
low-energy expression, eq. (M) . Indeed, this is one of the 
general problems to infer the full structure of xili^^n) 
entirely from the low-energy 'observable' x(g, w), and to 
settle ujo mayJpe a key point in the discrepancy between 
Radtke et alB and MP.B As it is not our purpose to 
discuss this point further, for the time being, we assume 
eq. (||) up to the cutoff energy uio in eq. (7), and arbi- 
trarily set uja = 0.4eV~ 4.6 x lO^K, following MP.I3) This 
value is used throughout in the following. Then we ob- 



the eigenequation (|T|), i.e., the vertex correction xf ^ tain .g^ = 0.57eV2 for Tc = 90K, stiU in disagreement 
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Fig. 3. The diagram (a) for the vertex corrected pairing potential 
yi^> and (b) for vi^'. 



with g2 ^ 0.41eV2 of MP.Q) This is not due to the size of 
a square lattice, as we see in Fig. 2, where Tc is shown as 
a function of .g^. At T = T^ = 90K, the 16x16 lattice is 
large enough for us to conclude g^ = 0.57eV^ in the self- 
consistent calculation including the effect of S(p, ia;„). 
A close inspection indicates that the disagreement orig- 
inates in details of xiQ^i^m)- In fact, we find g^ = 0.34, 
smaller than g'^ = 0.41 of MP, if we adopt the second line, 
instead of the third line, of eq. (7) for n = too. This 
means that Tc depends sensitively on how we prepare 
;\;(q',ii/„) in the low-energy regime. This is complemen- 
tary to the above remark on the high-energy contribution 
to Tc. As this quantitative difference is not of our primal 
concern either, deferring this problem, we choose to use 
our own definition, i.e., x(q, itj„) for n = is specified 
separately by eq. (7). The qualitative results presented 
below are not affected by this choice. 

Now let us discuss how we evaluate the vertex cor- 
rection. The pairing potentials Vv (p, iaj„;p', ia;„') and 
Vc (p, iw„;p',ia;n') including the vertex correction are 
diagrammatically represented by Fig. 3(a) and Fig. 3(b), 
respectively. These potentials at low frequencies uj„ = 
LOn' — t:T Ifpr n — n' = 0) are particularly studied by 
Monthoux.o' To see the effect on Tc precisely, however. 



-(2)/ 



we have to evaluate the kernel Kl ' (iu;„ , ia;„/ ) , eq. (|16D, 
for a full set of the Matsubara frequencies a;„ and ojn' , 
then the kernel must be diagonalized. In effect, this is 
not practical at present. Thus, as a tractable method, 
we set up perturbation theory to evaluate the vertex cor- 
rections to the eigenvalue k of the kernel. 

We shall make effective use of the results ob- 
tainable by means of the FFT. Let us introduce 
the eigenfunction ^^^^icun) for the largest eigen- 
value k'"'^-' of the kernel K'^^'>{iuJn,iuJn')', 'P''^\iuJn) = 
j:^,K(^')iiLu„,iu;n')^^^H^^n') = K^^^^^^KitOn) ■ Then the 
eigenvalue k including the vertex corrections is given by 



.•(1) 



c(2) 



.•(2) 



(20) 



where 



^(2) ^J2^^'H^^n)KPiiu^n,iu^n')'^^'H^^n'). (21) 

We assume <^'^^-'(ia;„) is normalized. On physical 



Fig. 4. Tc as a function of g^. Triangles; calculated without 

S(p, itiJn). Circles; including the effect of T:{p,iu!„). Diamonds; 

including the effect of S{p, iLO„) as well as the vertex corrections 

(2) (2) 

V„ and Vc . Two squares at T = 90K and 45K are calculated 

with Ti{p,iuJn) and V} but without Vv ■ 



grounds, the norm of ^'•^)(iw„) decreases quite rapidly 
as \iUn\ increases. Therefore, the sum over the Matsub- 
ara frequencies in eq. ( |2l| ) is allowed to be restricted 
in a narrow region around (n,n') ^ (0,0). In effect, 
we evaluate K^ (iu;„,iu;„/) for a 16 x 16 mesh around 
the Fermi energy. Moreover, in the remainder of the 
paper, the results are calculated on a 16 x 16 square 
lattice. Measured in terms of the weight |^^^-'(iw„)| , 

we find E|..„|<i5^T I^^^Hit^n)!^ = 0.98, 0.91 and 0.78 at 
T = Tc =90k745K and 22K, respectively. Even at low 
Tc, the error involved is not appreciable, for the coupling 
constant itself is small there. As Fig. 2 shows, a 16 x 16 
mesh in the momentum space is large enough to grasp 
the qualitative features caused by the vertex corrections. 
To prepare V^ (p, iwri;p',i'-^n') is most time- 
consuming. Therefore, first we use the bare Green's 
function G'^°^(p, iw„) instead of G(p, ia;„) to provide 

thus 
p 
calculated beforehand and G(p, iijj„) of the solution of 

eqs. (^ and (||), we calculate K^ ^(ia;„,iw„') in eq. (p^). 
Then, to evaluate k] , eq. (pi]), is straightforward, and 
the critical coupling g"^ at Tc is determined. Results thus 
obtained are shown in Fig. 4, where the triangles (with- 
out I](p, iw„)) and circles (with I](p, ia;„)) denote the 
results without the vertex corrections. (See Fig. 2). The 

(2) 

diamonds include the vertex correction Vv as well as 

(2) (2) 

Vc , while only the effect of Vc is taken into account 
for the two squares at T = 90K and 45K. 

Several points are noted from the figure. In the first 

(2) (2) 

place, both the effects of K and Vc are attractive 
on the whole, or enhances Tc fii the d-wave instability. 
The effect of Vc (Fig. 3(b))Xj however, is negligibly 
small, as noted by Monthoux.cJ) On the other hand, the 
effect of Vv (Fig. 3(a)) is prominent. In particular, it 
affects the result of MP, denoted by the circles interpo- 
lated with the solid line in Fig. 4, that the maximum 
transition temperature attainable in this model is about 



y/^^(p,iw„;p',itj„0- With (v^.^^^(p, itj„;p',icj„/) 
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Fig. 5. At T = 90K, k'-"^' and k'^) + k^^^ are shown as a function 
of g^. Squares, including only the effect of V} , are overlapping 
with circles to denote k'^' without the vertex corrections. 



lOOK.o' In fact, Tc as a function of g^ shows no sign 
of saturation, and keeps increasing beyond 200K when 

(2) 

the vertex correction Vv is taken into account. In this 
regard, the vertex correction has an effect more than a 
mere scale- up of the effective couphng constant g^ . 

Next, the effect of E(p, iw„) on V^ has to be inves- 
tigated. To this end, G(p, ia;„) to meet eqs. (||) and (g) 
is used to evaluate V^ (p, iw„;p',ici;„/). The maximum 
eigenvalues calculated for T = 90K arc shown in Fig. 5 
as a function of g^. The effect of I](p, ia;„) is to weaken 
the vertex corrections. The effect, however, is not ap- 
preciable for Tc — 90K, as we see from Fig. 5 in which 
we see g^ — 0.36 while we have g^ — 0.32 in Fig. 4 
in the case including the vertex corrections. Compar- 
ing these with g^ = 0.57 without the vertex corrections, 
we conclude that the correction due to S(p, iujn) in V^ 
is not important at least at Tc — 90K. In other words, 
if the coupling constant g should be evaluated to ac- 
count for Tc, our result is that the vertex. correction is 
not negligibly small ajL tlus temperatureJHS' at variance 
with previous results. cJ'B The discrepancy may be due 
to a high-energy contribution included in our calcula- 
tion, or it is traced back to the above finding of a slight 

(2) 

renormalization effect on V^ . On the other hand, for 
Tc = 180K, we find that g^ = 0.73 in Fig. 4 is modified to 
g^ = 1.57. The large modification in this case is due to 
a large coupling constant to realize that high transition 
temperature. The results in this regime must be taken 
with care. 

To the extent that the vertex corrections that we 
found for the pairing potentials are not negligible, the 
vertex corrections to eq. (|l|) should have to be investi- 
gated next.CJ-' The latter effect on I](p, iw„) will reduce 
Tc somewhat particularly through the pair propagator 



\G{p' ,iujn')\ in eq. (|16|), according to the above note, 
as a result of enhanced quasiparticle damping. There- 
fore, we will be ultimately led to a convergent result of 
Tc{g'^), somewhere in between the dashed and solid lines 
of Fig. 4. The results, however, would then indicate that 
Tc ~ lOOK is on the verge of practical applicability of 
this kind of perturbation theory in g^, as inferred from 



Fig. 4. Note that, for us in this context, to suffer a small 
correction is more important than to find out a high Tc- 

In summary, a result of this paper is presented in 
Fig. 4, though the result at high temperature is some- 
what modified as stated above. Applying perturbation 
theory to the eigenvalue of the kernel KiiujnjiuJn'), we 
estimated the vertex corrections to Tc as a function 
of the coupling constant g^ on the basis of the spin- 
fluctuation model of the high-Tc superconductivity. We 
found that the effect of Fig. 3(b) is numerically negligi- 
ble as far as the d^'^-y'^ pairing instability is concerned, 
while Fig. 3(a) enhances Tc appreciably. For Tc ^ lOOK, 
the effect of E(p, ia;„) mainly comes in through the pair 
propagator \G{p' ,iuJn>)\^ , dressing the vertex functions 
is not so important. In a strong-coupling regime at 
high temperatures, the vertex corrections become even 
qualitatively important, particularly in case where Tc in 
the one-loop Eliashberg calculation is substantially sup- 
pressed by lifetime effects. 
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